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t^j- | We investigate numerically and analytically the existence of thresholds for the pro- 
duction of Sine-Gordon single and multiple breathers from simple initial pulses. 
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1 Introduction 



The Sine-Gordon equation arises in physical problems such as one dimensional dislocations 
[1], long Josephson junctions[2], as well as in the mathematical treatment of constant negative 
curvature metric spaces [3]. 

Solitary wave solutions of the Sine-Gordon equation carry a topological winding number 
q. The solitary waves conserve their properties upon collision and are therefore called topo- 
logical solitons. The winding number zero sector q=0 consists of bound soliton-antisoliton 
solutions, the breathers, unbound soliton-antisoliton pairs[l], and phonons. The latter are 
small amplitude solutions, corresponding to the Klein-Gordon linear limit of the equation. 
The q=l sector solitary wave is the shelf or kink soliton. 

The Sine-Gordon equation is one of a gallery of nonlinear equations amenable to treat- 
ment by the Inverse Scattering Transform (1ST) method. The 1ST produces soliton solutions 
based on scattering data of Dirac-like equations for a potential connected to the sought solu- 
tion. The wave functions and the scattering data determine the potential - hence the name 
inverse scattering transform - through the integral equations named after Gelfand-Levitan 
and Marchenko [4] . 

The power of the 1ST method lies in the mapping of a nonlinear problem to an eigenvalue 
Sturm-Liouville problem. We recently used the 1ST method to treat the decay of distorted 
kinks to wobble states and phonons. [5] In the present work we investigate the decay of a 
square pulse to single and multiple breathers and phonons. The 1ST technique determines the 
constraints on the parameters of an initial pulse for it to decay into breathers and phonons. 
These constraints are satisfied at the location of a pole in the transmission amplitude of the 
1ST scattering problem. 

The decay of simple pulses to breathers was addressed some time ago by Kivshar and 
coworkers [6, 7]. Kivshar et al.[6, 7] studied the creation of breathers in long Josephson junc- 
tions and easy plane ferromagnets with an inhomogeneous Sine-Gordon equation. With a 
specific choice of the inhomogeneous source term and under the approximation of a negligible 
initial field configuration, the Zakharov-Shabat (ZS) equations[8, 9] are simplified consider- 
ably and an analytical solution of Kaup[10] could be implemented. In the present work we 
will focus on the homogeneous Sine-Gordon equation with a non negligible initial pulse. 

In section 2 we treat the decay of square pulse analytically using the AKNS[9] version of 
the Zakharov-Shabat [8] equations. The results are compared to a full numerical treatment 
of the equations with a more realistic pulse. Section 3 provides numerical support to the 
analytical results of section 2. Conclusions are drawn in section 4. 

2 Decay of a square pulse in the 1ST 

The Sine-Gordon equation is usually expressed in two ways. In space-time coordinates 
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and in light-cone coordinates x = -(X + T) t = -(X — T) with unit speed of light c = 1, 
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The Klein-Gordon equation for a particle of unit mass is obtained by linearizing Eqs.(l,2) 
Eq.(l) can be derived from the lagrangian 
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possesing a conserved energy 
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The 1ST method uses solutions of the Sine-Gordon and other nonlinear equations as a 
potential in a linear scattering problem of an ancillary wave function. The scattering problem 
is split into a time independent set of equations and a time evolution set. For the space-time 
formulation of the Sine-Gordon equation the equations are [9] 
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where £ is a complex spectral parameter. 

In eqs.(5),the X,T indices denote partial derivatives with respect to the physical coordinates. 
Consistency between the set of the first two equations and the second two equations in eq. (5) 
with a time independent £ yields the Sine-Gordon equation for u(x,t). The solitons of the 
Sine-Gordon equation correspond to the bound state sector of the spectrum of the AKNS 
equations [4]. These bound states appear as poles in the transmission coefficients. 

Two sets of independent solutions of the AKNS [9, 4] equations are defined. The first set 
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The second set ■?/> = | 1 | -0 = I Ul 

V ^2 / V ^2 

with boundary conditions at X — > oo 



e"^. (7) 
The two sets are connected by the scattering amplitudes a, a, b, b 

<j> = a(Z)$ + b(£)iP, 

= -a(0i> + b(0i>. (8) 

Poles in the transmission amplitude a(£) with Im(£) > generate the solitary waves. 
The problem of an initial square pulse can be solved analytically. The analytical solution 
can then be compared to the solution of eqs.(5) for more realistic pulses. 



Consider a static square pulse 
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u x = A (5(X)- S(w-X)), 



(9) 



where 6 denotes the Heaviside step function, w is the width of the pulse and A its height, 
and 6 represents the Dirac 5 function. The initial pulse is initially at its maximal (minimal) 
value. 

For 0, the solutions of eq.(5) are 



Vl = Q(- x ) e- iri{Xfl) + Q(x)Q(w - x)(B e 13 x + C e~ l 13 x ) + Q(x - w) Fe"* ', 
V2 = Q(x)Q{w-x){De tpx + Ee- ipx ) + Q{x-w) Ge iT,(x ' 0) , (10) 
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where B,C,D,E,F,G are £ dependent constant complex amplitudes and rj(X, 0) is defined in 
eq.(6). These amplitudes are determined by inserting eq.(10) in eq.(5) and using 0(0) = - 
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Poles in the transmission amplitude correspond to zeroes of the amplitude a in eq.(8), and 
to the vanishing of the amplitude P. The location of the poles is consequently determined 
by the implicit equation 



P = e^ w B (1 - 5 2 + 25 1 —-!-) + e~ l(3w C (1 - 5 2 + 25 ' J ^ ) = 0. (12) 

Single breather poles appear in pairs at £ = ±a + i(3, with (3 > to insure the convergence 
of the Neumann series of the Volterra integrals of the 1ST. [4], For a single breather at rest 
the spectral parameter obeys 

\£\ 2 = a 2 + (3 2 = 1 -. (13) 

An initially static pulse can decay to two breathers also. Momentum conservation requires 
the breathers to be emitted back to back with opposite velocity. A new parameter enters the 
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problem, the velocity of the breathers. For a breather in motion with velocity v the constraint 
of eq.(13) changes. It is easier to find the modified constraint in the light cone formulation. A 
Lorentz boost of a breather (a scalar field) amounts to a Lorentz transformation of the space- 
time arguments. In light-cone coordinates with unit speed of light, a boost with velocity v 
in the positive X direction transforms the light cone coordinates 




(14) 

Inspection of the spatial set of the Zakharov-Shabat equations in light-cone coordinates for 
the Sine-Gordon case[9] 



■ > 1 

Vl, x = -l&i - - U x V 2 , 

V2,x = i£v 2 + \u x v 1 , (15) 



shows that the spectral parameter for a moving breather transforms as £ — > \ £. The 

V 1 — v 

1 — v 

constraint of eq.(13) now reads a 2 + 3 2 = — -. As v — > 1, one pair of poles moves 

v ' 4 (l + v) 

towards the real £ axis while the other pair is pushed up to infinity. We here focus on the 
single breather production case and the double breather production is dealt with numeri- 
cally by integrating the Sine-Gordon equation(see section 3). The analytical treatment of the 
double breather case by means of the Hirota method [11] will be addressed in a future work. 
For the three breather production from a stationary pulse, one of the breathers remains at 
the location of the initial pulse. Therefore, three breather production can be treated as the 
single breather case. 

Eq.(12) was derived for the sharp edge square pulse of eq.(9). This pulse has infinite 
energy as may be seen by inserting eq.(12) in eq.(4), due to the 5 function discontinuities 

at X — 0, ±- w. A more realistic profile that still brings about the features of the square 

profile consists in the addition of lateral wings, continuous with the square pulse and its 
derivative at its edges. The modified initial pulse (centered at the origin for convenience) 
reads 



u(X,T = 0) = A {e~ z * 0(-i w-X) + e(X + ^ w) 0(^ w - X) + e~ z 'e(X - ^ w)), 
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with Sw, a width parameter for the wings. 

For eq.(16) we could not provide analytical solutions. Eqs.(5) were therefore solved nu- 
merically in order to verify the predictions of eq.(12). Eqs.(5) are firstly simplified by the 
substitution 



Ui -> v l e~ ip{x \ 
v 2 -> v 2 e lp{x) , 

P {x) = ^ x _ ^osHx)) d x y 

The spatial set of the AKNS equations(5) now becomes 
v, X = -![§£ + 

v %x = (f ( sin(u) + \[^ + §^]) v, e" 2 <**>, (18) 

with boundary conditions at X — > — oo, t> i — > 1, v 2 — > 0. 

We solved eqs.(18) using a fourth order Runge-Kutta method with double precision in 
order to achieve a maximal error of 1%. The spatial grid was fixed at 5X = .0002 with 
— 10 < X < 10. Poles in the transmission matrix were located by searching for D = 
abs(vi(-^> oo)) < 10~ 5 . For the wings width parameter we took Sw = 0.05, a small enough 
value to compare to the analytical results of the square pulse above. 

Figure 1 shows area plots for breathers production. The single breather regions were 
found using the analytical formula of eq.(12), curve labeled AKNS analytical, solutions of 
eq.(18), curve labeled AKNS numerical, and integration of the Sine-Gordon equation, curve 
labeled Sine-Gordon numerical. The thresholds for two breathers were obtained from the 
Sine-Gordon equation integration, while three breather thresholds followed from both the 
AKNS numerical treatment of eq.(18) and the square pulse formula of eq.(12). The curve 
labeled Sine-Gordon integration does not match the AKNS prediction. This is due to the 
fact that the observability of the breathers demands a minimal amplitude for them such 
that they do not blend into the phonon continuum. The thresholds found by the AKNS 
equations correspond to thresh — 0.5 + % 0, hence a vanishing breather height. The Sine- 
Gordon integration curve was drawn for breather amplitudes corresponding to a pole at 
£ = 0.499 + % 0.0316 at which the breather becomes clearly visible. We searched for poles in 
the parameter region < A < 10, 0.5 < w < 10, so that 5w stayed reasonably smaller than 
the width of the pulse. 
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Figure 1: Area plot for the parameter space leading to the production of Sine-Gordon 
breathers 
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Figure 2: v\ and v-i of eq.(18) as a function of distance for a value of £ correponding to a 
pole in the transmission matrix. 



The agreement between the square pulse predictions and the full AKNS treatment is 
satisfactory. From the plot we see that there is no single or triple breather production unless 
w >~ 7r. Figure 2 depicts a typical solution of eqs.(18) at a value of the spectral parameter 
corresponding to a pole. The wave functions are quite constant in each region as in the 
ansatz of eq.(10) corresponding to the square pulse. 

In the next section we depict the various breather production thresholds, and compare 
the breather profile to the theoretical expression. 

3 Decay of a square pulse, numerical simulations 

The numerical integration of the Sine-Gordon equation, was performed implementing the 
numerical code described in [5]. The spatial and temporal steps needed to achieve 1% 
accuracy in the final energy were dx = 0.02, dt = 0.01. The spatial extent of the integration 
region was taken to be X max = 600, with vanishing fields at infinity. In the examples depicted 
below we used 5w = 0.05. 

In section 2 we found two thresholds for single breather production from a square pulse. 
The width threshold reads w t hresh ~ tt- The height threshold can be obtained from the curve 
in figure 1. Figure 3 depicts a typical case below the width and height thresholds. A trail of 
phonons emanating from the origin is clearly observed, as well as the remnant of the original 
profile. The angular frequency of the oscillation of the central peak is greater than the upper 
bound for the breather pole a = ^ = 0.523 > 0.5, with T, the oscillation period. Moreover, 



9 




-90 -60 -30 30 60 90 



T=98.7 




-90 -60 -30 30 60 90 



Figure 3: Sine-Gordon solution for an initial pulse with A = 2.5 w = 2.5, below the width 
and height thresholds. 
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Figure 4: Sine-Gordon solution for an initial pulse with A = 1.5 w = 3.5, above the width 
threshold and below the height threshold. 
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the central peak does not have the shape of a breather. Figure 4 depicts a case below the 
height threshold and above the width threshold. Figure 5 shows a case above both width 
and height thresholds. In figures 3 and 4 only phonons are seen to emanate from the central 
region, while figure 5 shows a prominent peak centered at the origin and a small trail of 
phonons in both directions. From the graph of figure 5 the oscillation period is found to be 

I 7T 

T=6.6. Using £ = a + i/3, with (3 = y 0.25 — j3 2 , and a = — , we can obtain the location of 

the pole and evaluate eq.(12) with the amplitude and the height of the pulse as in figure 5. 
With these data, eq.(12) yields P = 2 1CT 3 , a value close to zero, as compared to P ~ 1 for 
£ far from a pole location. 

We can now compare the shape of the central peak in figure 5, presumed to be a breather, 
with the theoretical expression 



u(X,T) = 4 tan' 1 j, 
U = cos(2a (T + T )), 

V = cosh(2/3 X). (19) 

This correspondence is depicted in figure 6. Except for small phonon contributions, the 
profiles fit nicely. The central peak is a breather. 



11 



1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



a = 0.47656 (3 = 0.15126 




-25 



-15 



15 



25 



Figure 6: Central peak of the Sine-Gordon solution of figure 5 compared to eq.(19). 
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Figure 7: Sine-Gordon equation solution for an initial pulse with A = 4.1, w 
the single breather width threshold. 



2.5, below 



12 



T= 194.2 



-180 -130 



-30 



20 



70 



120 170 



Figure 8: Sine-Gordon equation solution for an initial pulse with A = 5.5, w = 3.1, above 
the single breather width threshold 

Figures 7 and 8 show numerical simulations for the two breather region of figure 1. A 
pair of back to back breathers receding from the center is observed. 
Figure 9 shows a three breather situation. 



4 Conclusions 

In the present investigation we considered the decay of an initially nonvanishing static square 
pulse. The results complement those of Kivshar et al.[6, 7]. However, contrarily to Kivshar 
et al. [7] we could not find any case for which a static square pulse leads to soliton-antisoliton 
pair production regardless of its width, height or energy. A different kind of initial pulse or 
a time dependent profile is apparently needed. 

Among the rich structure of the decaying pulse, the phonon spectrum certainly deserves 
attention. We found that in general the main features of the spectrum can be reproduced 
with the long time analytical expressions of Segur and Ablowitz[13], however, the fine details 
do not match existing expressions in the literature. 

We only touched upon the double breather problem and will deal with it in the future. 

We conclude with a remark concerning the experimental relevance of the present work. 
The initial pulse configuration we chose may be realized experimentally in long Josephson 
junctions of regular or high T c superconductors by injecting two oppositely directed bias 
currents at a distance w corresponding to the pulse width. By regulating the current intensity 
and the distance between injection points a well defined breather with a fixed oscillation 
period may be produced, with possible technological implications. 
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Figure 9: Sine-Gordon equation solution for an initial pulse with A = 8.2, w = 4 
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